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SUlitAWY 

This  paper  presents  a  simple  method  of  allowing  for  un~ 
certainties  in  thie  constant  tenns  (i.e.  right  hand  sides)  of 
a  linear  programming  problem^  and  hence  producing  realistic 
safety  margins  in  the  solution.  This  is  done  by  fitting  a 
mixture  of  uniform  distributions  to  the  assumed  distributions 
of  these  right  hand  sides,  and  using  a  particular  quadratic 
programming  algorithm.  C  ) 


M  --  Zt  1*  mil  knpan  tMt  llnMr  progrt— Ing  ha*  an  uneoapr^ 
■iting  addietion  to  booie  oolutlono.  This  noons  thot  it  will 
orto  uo  to  "go  oil  out"  for  our  objoctloop  ond  eonnot  bo  uood 
to  ootlaoto  tho  "oofoty  norgino"  thot  ofton  hovo  to  bo  protridod 
in  prootiool  oituotiono  to  guord  agoinot  uneortolntioo  of  ono 
sort  or  onothor.  go  OMitsig  (1956)  polntod  out*  otoehootie 
linoor  progri— ing  io  tho  thoorotieol  onouor  to  this  probloop 

I 

but  it  dooo  not  soon  to  hovo  boon  oppllod  oxtonoiooly.  This 
io  probably  bocouoo  tho  eonputotionf  Inoolvod  oro  uouolly  rothor 

hoofy*  Ofid  tho  ovoiloblo  data  on  chanet  offocto  aro  rartly 

D 

^Ptalaaaaaacb  to  ou^wrt  a  hoary  eonputation. 

^10  porpeoo  of  this  popor  io  to  show  that  an  laportant 
gppoial  olaoo  of  ouch  problono  can  bo  tolvod  fairly  coolly 

*a  *!•  ' 

using  a  otandard  quadratic  prograaning  algorlthn.  go  for  oo 
1  lOMOf  this  algorithn  haa  not  yot  boon  prograaotd  for  o  conputor 
but  1  hops  that  it  will  bo  in  tho  not  too  distant  future.  It 
Oliaald  tnon  bo  tiorth  uiiilo  to  food  in  cotloatoo  of  tho  vin- 


oortointioo*  and  thoir  oonooqponeoop  into  tho  cooiputotlono 
ovon  if  thooo  ootlMtoo  aro  quite  crude. 


fho  elaoa  of  problono  eonoldorod  here  has  boon  dlocuoood 
by  nany  authorot  Oanttlg  (1953) »  Beale  (1933)#  fbrguoon  ond 
Oantaig  (1936)*  Unaghraby  (1939)#  (1960)#  Oantsig  and  lOidonolcy 


'  1. 


2 


(1960)f  and  Nadansky  (19^0)> 

The  nonttochastic  version  of  the  model  contains  equations 
of  the  form 

(1)  six  -  bj. 

where  is  some  linear  function  of  the  variables  of  the 
problem,  denoted  collectively  by  the  vector  x,  and  where 
is  a  number  that  one  pretends  Is  known  exactly  when  it  really 
ia  not.  Vor  example,  it  may  be  the  volume  of  sales  in  sosie 
future  time-period.  The  problem  is  to  choose  a  nonnegative  ;|c 
satisfying  all  these  equations  so  as  to  minimise  some  linear 
function  c'x. 

Zn  the  stochastic  version  of  the  model,  the  b^  are  regarded 
as  random  variables  with  known  probability  distributions,  ^e 
vector  ^  has  to  be  chosen  before  we  know  the  actual  values  of 
the  b^;  but  we  are  no  longer  required  to  satisfy  the  equations 
(1)  exactly.  Instead  we  have  a  penalty  of  f^(2  0)  for  each 
unit  by  which  g^x  falls  short  of  b^,  and  a  penalty  of  f^(2  0) 
for  each  unit  by  which  g^x  exceeds  b^,  to  be  added  to  the  direct 
cost  c';|c.  This  penalty  may  represent  a  tangible  loss,  such  as 
failure  to  sell  one's  goods  because  the  demand  is  inadequate, 
or  an  intangible  loss  of  goodwill  caused  by  failure  to  meet  a 
requirement,  or  a  combination  of  the  two.  The  problem  is  then 
to  Choose  a  nonnegative  x  that  minimizes  the  expected  total  cost 


(2) 


where  denotes  the  neen  value  of  the  penalty  Incurred  on 
the  i-th  equation.  To  represent  rstheinatlcally,  we  write 
(1)  in  the  form 

(5)  *i«  -  •  ‘’i* 

Where 

(^)  ^i  •  y?  *  ^1*  ^1  2  ^  0*  ^ 

and 

(5)  Pj  -  I(fj  Pj  ♦  fj  yp. 

...  , 

IT;  In  practice  there  will  often  be  no  penalty  coat  if  y^  >  0* 
^  ao  that  «  0.  The  nonstochastic  model  would  then  normally 
appear  in  the  form  2 

Zn  fact  it  is  convenient  to  reduce  the  general  problem 
ta  one  with  f^  ■  0#  by  adding  f^  times  equation  (J)  to  equation 

(a)*  Ve  then  have 

(6)  T--rfJ  b74 

idiare  denotes  the  mean  value  of  b^»  and 

(7)  p*  -  l((fj  ♦  i^)yp. 

Zt  ia  perhaps  worth  noting  that,  although  we  have  assumed 
that  and  f^  are  both  nonnegative,  we  really  only  require 
'  that  f ^  ♦  f7  be  nonnegative  to  Justify  the  siathenatics. 
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Given  a  problem  with  ■  0*  it  is  convenient  to  write 

(8) 

and  to  refer  to  as  the  safety  margin,  or  margin,  in  the  i*th 
constraint.  The  essential  contribution  of  the  stochastic 
feature  of  the  model  is  to  provide  a  means  of  estimating  proper 
values  for  these  margins.  In  the  nonstochastic  model  they  are 
zero  on  all  the  operative  constraints. 

Another  general  point  to  be  noted  is  that  the  correlations 
between  different  b^  do  not  enter  the  problem.  They  have  no 
effect  on  the  mean  value  of  the  penalty  cost  —  although  of 
course  they  could  have  an  appreciable  effect  on  its  variance. 

3»  mSTIIIQ  WBTH0D3  OP  SOLOTIOW 

Methods  for  solving  these  problems  are  known.  If  all  the 
probability  distributions  are  discrete,  any  2-«tage  linear 
programming  problem  can  be  reduced  to  a  vast  linear  programming 
problem,  which  can  nevertheless  be  solved  using  the  decomposition 
principle  of  Dantzig  and  Wolfe  (i960),  as  Dantzig  and  Madansky 
(i960)  have  pointed  out.  The  simple  type  of  penalty  coat  con¬ 
sidered  here  is  easier  than  this.  It  leads  Immediately  to  a 
convex  separable  nonlinear  objective  function,  which  can  be 
handled  by  the  metnods  of  Chames  and  Lemke  (195^)* 

Rlmaghraby  (i960)  has  suggested  that  a  continuous  dis¬ 
tribution  function  could  be  handled  using  the  Lagrangian  differ¬ 
ential  gradient  method  of  Hurwioz  (1957)*  But  the 


Halted  cooputatlonal  •jqptriMce  with  th«te  Mthods  Is  not 
voxy  soeourocinc- 


His  ttartiiis  point  for  tho  approach  proposed  hare  la  the 
fact«  noted  bgr  Oaatzlg  (1953)*  that  If  the  randoa  variable 
has  a  uniform  (rectangular)  distribution  over  soms  specified 
range  then  the  objective  l\uiction  is  quadratic.  This  follows 
from  the  fact  that  (assuming  f^  ■  0)»  d?^/dy^  equals  minus 
times  the  probability  that  <  0*  in  fact  we  find  that  if 
b^  is  unifonUy  distributed  between  —  r.  and  B.  •¥  then 


P^.O 


Pi  1  «“l. 


L*  5rJ 


♦sir" 


1  <  Pi  <  ‘•i 


Pi  1 
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This  function,  plotted  against  y^,  look*  as  follows. 


And  It  turns  out  that  this  function  can  be  represented  very 
simply,  by  wrltlrig 


(9)  -  '•i  ♦  *11  -  *21  -  *51. 

Xu  >  0,  >  0,  Xjj^  >  0, 

(10)  -  r“  (x,^j  + 

In  other  words  the  constraints  (1)  or  (3)  are  replaced  by 


(11) 


ilx  -  Xu  ♦  Xu  +  Xj,  -  Ej  ♦  rj, 


and  (2)  Is  replaced  by 
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(12)  *  •  g')"  ♦  >^l/‘*ri)- 

To  verify  theie  formulae,  first  consider  the  situation 
when  d****  that  we  can  put  ■  0* 

and  •  0,  as  it  should  be.  Now  suppose  that  y^  <  r^.  Then 
we  cannot  put  x^^  ■  ■  0,  and  F^  will  have  to  be  greater 

than  zero.  But  we  would  be  crazy  not  to  put  x^^  ■  0,  since  t 
positive  x^i^  would  involve  an  unnecessarily  large  x^^  or 
at  SQise  extra  expense  and  no  benefit.  So  we  put  ■  0  and 
make  up  the  difference  between  y^^  and  r^^  by  x^^  -f  How 

iF^/ixg^  ■  f^,  and  •  ^l*31^*'l’  ** 

increase  x^^  <e  it  is  cheaper  to  increase  x^^  if  x^^  <  2r^, 
^  and  otherwise  it  is  cheaper  to  increase  Hence,  as  y^ 

decreases  from  r^j  F^  increases  quadratically  to  when 
•  -r^.  and  thereafter  f^  increases  linearly.  In  fact  f^ 
represented  faithfully  in  all  cases. 

If  it  is  more  convenient,  we  can  handle  a  problem  with 
B;  ^  ^  0  directly,  by  writing 

W  (W)  -  fjxij  ♦  <7*21  -  fj»5i  * 

■ ' 

The  problem  of  minimizing  T  subject  to  these  constraints 

[r.  is  a  typical  quadratic  programming  problem  in  which  the  number 

t 

\  of  quadratic  terms  equals  the  number  of  constraints  with  random 
right  hand  sides,  which  is  almost  certainly  much  smaller  than 

I 

the  total  number  of  variables  in  the  problem.  The  problem  is 

k 


*  i 


therefore  suited  to  the  algorithm  outlined  on  pp.  235-536  of 
Beale  (1959)‘  This  Is  essentially  the  algorithm  presented  by 
Beale  (1955)  and  more  fully  In  Beale  (1959) •  only  dlfferenco 

Is  that  the  objective  function  Is  stored,  not  as  a  vast  square 
matrix,  but  as  a  set  of  (r<fl)  linear  forma  X^, 

It  being  understood  that 

T  ■  Xq  ••• 

So  we  can  deal  with  unlfomly  distributed  right  hand  sides 
fairly  easily.  In  some  situations  we  may  know  so  little  about 
the  true  distributions  of  these  right  hand  sides  that  such  an 
assumption  Is  as  good  as  any.  But  In  other  situations  this 
will  not  be  very  satisfactory.  For  example,  suppose  we  have 
a  production  scheduling  problem  In  which  sales  are  treated  as 
random  variables,  but  we  wish  to  lay  down  the  production  schedule 
In  advance.  T^ien  the  nonstochastic  model  will  contain  some 
constraints  of  the  form 


>  3^, 


Iq  ♦  ^1  ♦  >  Si  >  $2 


representing  tne  facts  that  Initial  Inventory  plus  total  pro¬ 
duction  up  to  the  end  of  any  time-period  must  not  be  less  than 
total  sales  up  to  the  end  of  the  time-period.  Now  we  might 
possibly  assume  that  has  a  uniform  distribution.  But  it 
would  hardly  be  consistent  to  assume  t^iat  ♦  Sg  also  had  a 


uniform  distribution. 
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A  posalbly  moi’e  fundlamental  objection  to  tne  uniform 
llstr-lbutlori  la  tiiat  It  aasunea  tf.at  there  la  no  dar.^jer  t:At 
the  random  variable  will  differ  from  Its  mean  by  more  than 
1.7?2  standard  deviations.  And  yet  it  la  clear  that  In  some 
situations  wider  safety  margins  are  needed. 

But  now  we  observe  that  the  analysis  can  be  carried  through 
If  the  distribution  of  Is  a  mixture  of  a  small  number  of 
uniform  distributions  almost  as  easily  as  If  It  Is  a  single 
uniform  distribution.  Suppose  that  there  Is  a  probability  p^j 
that  b^  Is  unlforenly  distributed  between  ai* . 

•f  rj^j,  for  J  ■  1,  ...,  Then  If  we  replace  the  1-th  con¬ 

straint  by  the  constraints 

(i*)  ^5  -  +  Xjjj  ♦  Xjjj  .  ♦  r^j, 

we  have 

(15)  T  .  c'x  +  I  Ifj  Pjj  (Xj^j  ♦  *31/“ '•ij)- 

This  appioacn  would  be  unattractive  If  we  nad  to  use  a 
large  value  of  k^,  l.e.  a  laige  number  of  unlfonr.ly  distributed 
components,  for  each  constraint  with  a  random  right  hand  side. 
One  might  think,  for  example  that  several  components  would  be 
needed  to  approximate  a  nonnal  frequency  function  by  a  step 
function.  But  the  penalty  cost  Is  not  represented  by  the 
frequency  function,  but  by  the  Integral  of  the  cumulative  dis¬ 
tribution  function.  And  t.^te  process  of  Integrating  twice 
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smootns  out-  the  comers  In  the  approximation  In  a  most  gratifying 
way.  An  approximation  (fitted  by  the  method  of  momenta)  to 
the  nonnal  distribution  by  Just  two  uniform  distributions  gives 
an  expected  penalty  cost  that  Is  correct  to  within  1^  for 
any  margin  less  than  1.5  standard  deviations.  For  a  margin 
between  1.5  and  2.5  standard  deviations  the  approximation 
exaggerates  the  expected  penalty  cost  by  up  to  50!ll.  This  may 
seem  an  undesirably  large  exaggeration,  but  It  would  arise 
equally  If  one  used  a  genuine  normal  distribution  with  a 
standard  deviation  overestimated  by  5jK.  So  In  practice  It 
will  often  be  unrealistic  to  try  for  greater  precision  than 
trAt  given  by  two  uniformly  distributed  components.  Tf.e  de¬ 
tails  of  the  method  of  fitting  are  discussed  In  the  appendices, 
tne  results  for  a  normal  distribution  being  Illustrated  In 
Figs.  1  and  2. 

We  therefore  need  Just  one  extra  constraint,  and  two 
quadratic  variables  In  the  objective  function,  for  each  con¬ 
straint  In  tne  original  problem  with  a  random  right  hand  side. 

Since  this  will  usually  apply  to  not  more  than  half  the  constraints, 
this  Is  not  a  large  addition. 

5.  THE  RSLATIOWSHIP  WITH  CHAMC^OWSTKAIMKD  UMIAR  FFOOWAIWIIIO 

It  Is  of  Interest  to  compare  tnls  "2— stage"  approach  to 
the  problem  of  unknown  right  hand  sides  with  chance— constrained 
programming.  Introduced  by  C..ames,  Cooper  and  Symonda  (1958). 

Of  course  trils  problem  Is  not  tne  most  general  possible  problem 
for  either  ci  ance— constrained  or  2— stage  linear  programming; 
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but  it  Is  often  illuminating  to  compare  different  approaches 
with  reference  to  a  specific  simple  class  of  problems. 

One  difference  is  that  the  2-stage  program  minimizes  the 
total  costj  consisting  of  a  direct  cost  and  a  penalty  cost 
representing  the  average  level  of  failure  to  satisfy  the  original 
inequality  constraints.  In  chance-constrained  programming  one 
pre-selects  a  tolerable  long-term  failure  level,  and  simply 
miniaizes  the  direct  cost  without  exceeding  this  failure  level 
on  the  average.  Such  an  appi'oach  may  be  advantageous  if  the 
penalty  costs  are  hard  to  assess  quantitatively.  For  our  pro¬ 
blem  it  can  easily  be  incorporated  into  the  2— stage  model  by 
using  parametric  programming.  As  Beale  (1939)  Indicates, 
parametric  quadratic  programming  is  straightforward  as  long  as 
the  parameter  is  confined  to  the  right  hand  sides  and  the  linear 
part  of  the  objective  function.  To  apply  it  in  this  context 
we  must  therefore  keep  the  penalty  cost  function  fixed,  and 
gradually  scale  down  the  direct  cost  to  represent  a  gradually 
Increasing  relative  importance  to  be  assigned  to  failure. 
Corresponding  to  each  parameter  value,  one  can  coiTipute  the 
direct  cost  and  the  average  level  of  infeasibility,  measured 
by  the  average  penalty  cost  incurred.  One  can  then  select  that 
parameter  value  corresponding  to  the  selected  tolerance  for 
inftasibillty . 

But  the  more  fundamental  difference  between  2-6tage  and 
chance-constrained  programming  is  that,  for  any  given  values 
of  the  variables  and  the  random  elements,  the  infeaslblllty  is 
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measured  In  cnance-cons trained  programming  essentially  by  the 
number  of  constraints  that  are  violated,  and  not  by  the  extent 
to  which  they  are  violated.  (Actually  this  Is  an  oversimpli¬ 
fication.  Violations  of  some  constraints  may  have  to  be  re¬ 
garded  as  mor'e  serious  than  violations  of  others,  to  allow  the 
overall  probabilities  of  violating  different  constraints  to 
balance  In  an  optimum  program.)  Tne  question  of  which  of  these 
Is  the  more  reasonable  measure  must  depend  on  the  application. 
Madansky  (I960)  suggests  an  alternative  form  of  chance-con¬ 
strained  programming  In  which  the  Infeaslblllty  is  measured  by 
whether  or  not  all  the  constraints  are  satisfied  simultaneously. 
To  apply  this  to  our  problem,  we  would  need  to  Icnow  the  Joint 
distribution  of  all  tne  random  variables,  whereas  the  marginal 
distributions  suffice  for  the  other  approaches. 

It  Is  possible  to  compromise  between  our  formulation  of 
the  problem  and  cnance— constt'alned  programming,  by  accepting 
our  measure  of  Infeaslblllty,  computing  the  parametric  family 
of  optimum  solutions  for  different  relative  weights  on  cost  and 
feasibility,  and  then  selecting  from  this  family  the  solution 
with  a  specified  average  total  number  of  violated  constraints. 

In  some  circumstances  this  might  be  a  very  practical  procedure, 
vlvlng  an  easily  comprehensible  If  Indirect  control  over  the 
level  of  Infeaslblllty. 


kfnmix  1 

Afrnoxinknm  a  symiitkical  DisnuRfrioii  by  a 
Mrnms  of  otiifowi  nxaTWBonoiis 

To  apply  the  nethoda  of  Section  4,  the  aasumed  distribution 
for  the  random  variable  b  muat  be  approxlnated  by  a  mixture  of 
a  small  number  of  uniform  distributions.  It  Is  natural  to  use 
the  method  of  moments  to  define  this  approximation.  This  Is 
a  sosMwnat  arbitrary  aeclslon*  bvt  It  seems  natural,  because 

(a)  the  moments  of  the  given  distribution  can  usually  be 
readily  evaluated,  and 

(b)  It  tends  to  emphasize  the  tails  of  the  distribution, 
where  a  good  fit  Is  most  Important. 

The  problem  Is  greatly  simplified  If  b  has  a  symmetrical 
distribution.  It  Is  then  clear  that  we  should  make  t  e  mean 
of  each  component  equal  to  the  mean  of  the  distribution,  and  It 
only  remains  to  choose  the  proportions  p^  and  half-ranges  r^ 
to  be  allotted  to  the  different  components.  With  k  components, 
we  caa  fit  up  to  the  (4k  —  2)— th  aoment,  and  the  equations  to  be 
oatiefiod  ore  as  follows: 


see 

-  2) 


•  1 


3u 


2 


•  ~  ^^^4k  —  2* 
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where  dcnotec  the  1— th  ooiaent  about  the  mean  of  the  (syiaaetric) 
distribution  being  fitted. 

Writing  (21  ♦  “  '^l* 

equations  reduce  to 


.2k-l 


There  are  the  standard  formilac  for  the  weights  and 
ordinates  for  aauss-type  quadrature  forciulae.  For  convenience 
we  record  the  solutions  for  1:  •  1,  2  or  3. 

If  k  .  1,  -  1. 

If  k  ■  2,  X^  and  X^  arc  the  roots  of  the  quadratic  equation 


X  1 

1  .0, 

and 

Pi  -  (v^  -  Xg)/(Xj  -  X^) 

^2  *  ^'^1  ”* 

If  k  •  3,  X^,  X^  and  X^  are  the  x^ots  of  the  cubic  equation 


M4(J| 
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and 

Pj  -  (*8  -  Vj(X2  ♦  X3)  4  XgX3)/{Xj  -  Xg)(Xj  -  X3) 
P2  •  (*2  -  *1(^1  ♦  Xj)  ♦  ^i»3)/(Xj  -  XjXXj  -  X3; 

Pa  ■  <*2  -  \ih  +  ^3)  ♦  V2)/(^3  -  -  ^s). 


Pot  a  nomBl  dlatrlbutlon  with  variance 

»2  -  1.3o*,  V3  .  1.3.5.706^  .  1.3.5.7.90®, 

•nd  we  have  the  follcwln*  reaulto. 


Vj  •  1.3o®, 

•  1.3.5.7.9.110 


0.1838 

0.8162 

0.0154 

0.3446 

0.6400 


1.73210 

2.85700 

1.35560 

3.7504o 

2.3^ 

1.15440 


*»•  fitted  frwmency  nmctloos  are  plotted  in  W*.  1, 

•ad  the  reauitins  •vei«*e  ponalty  coata,  0  .  i,  f  . 

fl  -  0,  are  plotted  on  a  logarlthnlc  scale  in  Fig.  2.  it  win 
»•  aMa  that  the  aingle  unlfora  diatrlhutloo  (k  ■  1)  givea  a 
»aiy  ado«»te  fit  for  aaisma  less  thvi  o,  but  ««doreatli*tea 


F-2404 
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the  penalty  for  larger  margins.  The  mixture  of  2  uniform 
distributions  (k  »  2)  gives  an  adequate  fit  for  margins  less 
than  2.5o,  and  this  should  be  satisfactory  for  most  practical 
purposes.  The  Improvement  obtained  by  taking  a  third  component 
(k  ■  3)  l8  not  spectacular.  One  gets  an  adequate  fit  right 
out  to  but  there  Is  an  awkward  trough  In  the  curve  around 

2.3®  where  the  penalty  Is  underestimated  by  3^1^*  This  Is  per¬ 
haps  because  the  tenth  moment,  which  Is  used  to  fit  tnls  mix¬ 
ture,  gives  too  much  weight  to  the  extreme  ta^ls  of  the  dis¬ 
tribution.  A  better  fit  over  a  somewhat  shorter  range  of 
margins  could  perhaps  be  obtained  by  fitting  the  moments  of 
a  truncated  nonnal  distribution. 
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APfgWDIX  2 

APPROXINATIIK)  AN  UN!;YlttBTRICAL  DISTNIBOTION 

BY  A  mrnmi  of  mtirom  DisnimmoKs 

Tht  problem  of  fitting  •  mixture  of  unlfom  distributions 
by  the  method  of  moments  to  sn  unsymmetrlesl  distribution  seems 
to  be  much  more  swlcwftrd  then  fitting  s  sysmetrlcsl  distribution. 
We  now  have  to  choose  3  parameters  for  each  cosiponent:  the 
proportion  the  mean  and  the  half-range  r^.  With  k  com¬ 

ponents  we  can  therefore  fit  up  to  the  (3k  -  l)-th  moment. 
Without  loss  of  generality  we  can  assume  that  the  mean  of  the 
distribution  to  be  fitted  Is  zero. 

The  r-th  moment  of  the  mixture  Is  given  by 


Z 

1 


?rrtl7r^ 


((m^  ♦ 


). 


Xf  u  denotes  the  ]>>th  moment  about  the  mean  of  the  distribution 
r 

to  be  fitted,  we  therefore  obtain  the  following  equations  for 
fitting  2  components t 


e 


There  does  not  seem  to  be  any  neat  vray  of  solving  these 

equations.  But  we  can  tackle  them  by  solving  the  first  2 

equations  for  and  in  terrr.s  of  and  substituting 

these  values  Into  the  3-rd  and  4— th  equation  and  solving  for 
2  2 

and  r^  In  term?  of  and  m^,  and  then  substituting  all 
these  values  Into  the  5-th  and  6-th  equations  to  detennlne 
and  m2*  V/e  find  that 


(A2.2)  -  njg)*  Pg  -  -  1112)  V 

Jj 

(A2.3)  ^  3^2  ♦  2mj^m2  “  ”l^'  **2  "  tn^  ^2 

•J 

When  we  substitute  these  values  into  the  last  2  equations 

■4 

of  (A2.1)  It  la  convenient  to  write 


(A2.4) 


The  equations  then  reduce  to 

(A2.5)  ^10  ^llh  ♦  ^12^1 

(A2.o)  ^20  ^  ^21^1  ^  ^22^1  ^  ^23^1  * 

AlO  •  ^3  ^^2  -  -  12^2^  - 

All  *  ^3^ 

Ai2  •  ♦  ^^2* 


where 


*20  "  ♦  ®S*2 

Agi  -  -  5l^  -  12u2)|  -  16>| 

*22  ■  ♦  ***3*2 
*23  “  *  **!• 

W»  can  manipulate  these  equations  to  solve  for  In 
terns  of  X^,  and  obtain  the  fonmila 

(«.7)  .  «i  *  -  9»^K>*2  ^  (3U5  -  ^ 

♦  (5|4|j  -  9ii2)\) 

Finally,  substituting  this  back  Into  (A2.5)  ve  obUln  a 
sextlc  equation  for  Xg,  which  can  be  written  as 

(A2.8)  ♦  *2^  ♦  *3*2  ♦  >4*2  ♦  “s*!  ♦  '•6*l  ’  <>- 

Where 

•o  • 

*2  •  I6ai*||»5  -  25l^|t2  ♦  9ai|l*®|»*  -  76842u5  -  8luju^. 

*3  •  2*C»i5H*i»5  -  -  1254J  ♦  675l4uf  -  laOOijulii* 

-  121Si4u4  -  12811*  2160|i||i|  ♦  729|i|, 

•4  -  56»»|  -  l*»Ji2HjU5  -  JOOUguJ  -  32<»i|u^  ♦  108()»4u4 
720i|ii|  -  972n|, 


i>-e«d» 

Si'O 


.  igzujtij  -  ‘•(xv®  ♦  utonlnj,  -  -  129614, 

“6  "  256*4- 

The  suggested  computational  procedure  is  therefore  as 
follows  I 

1.  Compute  the  coefficients  of  (A2.8)«  and  find  a  negative  real 
root  of  the  equation.  If  there  are  no  negative  r'iots,  then  the 
problem  is  insoluble.  If  there  are  any,  then  there  must  be  an 
even  number.  It  seems  likely  that  the  smallest  in  absolute 
value  will  be  appropriate,  but  this  is  not  exactly  clear. 

2.  Given  the  value  of  deduce  from  (A2.7)* 

3.  Find  the  roots  of  the  quadratic  equation 

2 

X  -  X^x  4-  Xg  •  0. 

Tnese  are  m^  and  m2 • 

2  2 

4.  Deduce  p^,  P2,  r^^  and  from  (A2.2)  and  (A2.3)*  Since 

Xg  <  0,  mj  and  irig  must  have  opposite  signs,  so  pj  and  Pg  are 

2  2 

certain  to  be  positive.  But  r^  or  might  prove  to  be  negative 
indicating  that  the  problem  has  no  solution,  at  least  for  this 
value  of  X2 • 


w 


Average  p«noiry  cost 


Fig  2 — AppfOftimatiog  o  normol  distribution 
by  0  ffiixturt  of  uniform  distributions: 
Avtrogt  ptnolty  costs  f^ ■!,  ft  *0) 
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